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Abstract: Static excited-state polarisa- 
bilities and hyperpolarisabilities of 
three Ru" ammine complexes are com- 
puted at the density functional theory 
(DFT) and several correlated ab initio 
levels. Most accurate modelling of the 
low energy electronic absorption spec- 
trum is obtained with the hybrid func- 
tional B3LYP, B3P86 or M06 for 
the complex [Ru u (NH 3 ) 5 (MeQ + )] 3+ 
(MeQ + = A r -methyl-4,4'-bipyridinium, 
3) in acetonitrile. The match with ex- 
perimental data is less good for [Ru 11 - 
(NH 3 ) 5 L] 3 + (L = N-methylpyrazinium, 
2 ; A f -methyl-4-{£,£-4-(4-pyridyl)buta- 
1,3-dienyljpyridinium, 4). These calcu- 
lations confirm that the first dipole- 
allowed excited state (FDAES) has 
metal-to-ligand charge-transfer (MLCT) 
character. Both the solution and gas- 



phase results obtained for 3 by using 
B3LYP, B3P86 or M06 are very similar 
to those from restricted active-space 
SCF second-order perturbation theory 
(RASPT2) with a very large basis set 
and large active space. However, the 
time-dependent DFT 2 max predictions 
from the long-range corrected func- 
tional CAM-B3LYP, LC-wPBE and 
wB97XB and also the fully ab initio 
resolution of identity approximate cou- 
pled-cluster method (gas-phase only) 
are less accurate for all three com- 
plexes. The ground state (GS) two- 
state approximation first hyperpolarisa- 
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bility /3 2SA for 3 from RASPT2 is very 
close to that derived experimentally via 
hyper-Rayleigh scattering, whereas the 
corresponding DFT-based values are 
considerably larger. The /3 responses 
calculated by using B3LYP, B3P86 or 
M06 increase markedly as the it-conju- 
gation extends on moving along the 
series 2^4, for both the GS and 
FDAES species. All three functionals 
predict substantial FDAES /? enhance- 
ments for each complex, increasing 
with the it-conjugation, up to about 
sevenfold for 4. Also, the computed 
second hyperpolarisabilities y generally 
increase in the FDAES, but the results 
vary between the different functionals. 



Introduction 

The promise of diverse applications including optical data 
processing and biological imaging has stimulated much in- 
terest in organic nonlinear optical (NLO) materials.' 11 
Within this field, organotransition metal complexes offer in- 
triguing possibilities for creating new multifunctional materi- 
als in which potentially useful optical behaviour is combined 
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with the redox, magnetic and other properties characteristic 
of such compounds.' 2 ' As a means to enhance the prospects 
for molecular materials, approaches to modulating reversi- 
bly molecular NLO properties have attracted considerable 
attention recently.' 3 ' The first report of a very pronounced 
and reversible redox-switching of the first hyperpolarisabili- 
ty (the origin of quadratic molecular NLO effects) clearly 
demonstrated the potential significance of metal complexes 
in this area.' 4 ' This work has been extended to second har- 
monic generation (SHG) in Langmuir-Blodgett (LB) thin 
films,'" 1 ' and various related solution and molecular-level the- 
oretical studies involving both quadratic and cubic NLO ef- 
fects have been described.' 6 ' 

Although the switching of NLO responses via redox 
chemistry has attracted much attention, using other stimuli 
is also of great interest. Light-induced (photochromic) mo- 
lecular rearrangements have been explored relatively 
widely.' 3,7 ' However, the speed of switching accompanying 
such structural changes is often quite limited, and their ef- 
fects on macroscopic structures may be substantial. Populat- 
ing transient electronic excited states should allow much 
faster modulation effects without significant structural 
changes. Early CNDO/S calculations on some simple dipolar 
molecules, such as 4-nitroaniline, indicated that /3 responses 
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can be increased in excited states,' 8 ' and these predictions 
were verified subsequently.' 9 ' Various experimental studies 
with purely organic chromophores reveal similar behaviour, 
mostly focusing on second hyperpolarisabilities y,' 10 ' and at- 
tempts to determine excited-state Rvalues by hyper-Ray- 
leigh scattering (HRS) have been reported.' 11 ' However, 
such measurements are fraught with complications and 
hence there is currently no reliable experimental method by 
which such molecular-level excited-state responses can be 
determined. Therefore, the establishment of accurate theo- 
retical approaches is of significant interest, as a guide to em- 
pirical studies that will most likely focus on the switching of 
bulk effects like SHG. 

Some time ago, Sakaguchi et al. noted a small photo- 
induced switching of SHG at 295 nm from alternating 
and highly diluted LB films containing the amide-substituted 
[Ru"(2,2'-bpy) 3 ] 2+ (bpy = bipyridyl) derivative 1 (Figure l).' 12 ' 
This observation was attributed to changes in /? on MLCT 
excitation. However, the ground-state (GS) complex shows 
an intense 2,2'-bpy-based absorption near the SHG wave- 
length, so excitation-induced changes in absorption may 
affect the SHG signal. Notably, very few other related ex- 
perimental or theoretical studies with metal complexes have 
been reported to our knowledge (and these concern only 
y responses).' 13 ' 

The [Ru n (NH 3 ) 5 ] 2+ complexes that we have studied as 
redox-switchable NLO chromophores' 4,5 ' show intense 
MLCT bands in the visible region. Consequently, they show 
very large /S responses that compare favourably with those 
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Figure 1. A ruthenium complex (1) studied previously for ps timescale 
photoswitching of SHG in LB films,' 12 ' and the complexes 2-4 considered 
in the present work. 



of purely organic chromophores with moderate it-conjuga- 
tion lengths.' 2 ' 14 ' These systems are therefore highly attrac- 
tive subjects for photoswitching studies, and we describe 
here theoretical investigations, which indicate that MLCT 
excitations lead to large changes in molecular NLO respons- 
es. The use of time-dependent density functional theory 
(TD-DFT) and ab initio methods to predict hyperpolarisa- 
bilities of GS chromophores is now well developed,' 15 ' but 
relevant considerations of electronic excited states are re- 
stricted to polarisabilities and/or small molecules.' 16 ' 



Results and Discussion 

Our calculations focus on the series of complexes 2-4 
(Figure 1). Complex 2 is well studied,' 17 ' and Stark spectros- 
copy in a 1:1 glycerol/water glass has been used to derive 
a modest static first hyperpolarisability /3 0 (< 1(T 29 esu) for 
the salt [2][BF 4 ] 3 .' 18 ' We have used HRS in acetonitrile 
(MeCN) solutions and Stark measurements in butyronitrile 
(PrCN) glasses to afford substantially larger /3 0 values (HRS/ 
Stark, lCT 3(1 esu) of 246/240 for [3] [PF 6 ] 3 ' 19,201 and 744/1092 
for [4][PF 6 ] 3 .' 21 ' Note that we consistently use the so-called 
T-convention in this work,' 22 ' whereas previous reports' 19 " 21 ' 
use the B-convention; to convert to the T-convention, we 
use J3 T =2/3 B . In addition, all of our values refer to the domi- 
nant component along the long molecular axis, approximate- 
ly parallel to the dipole moment. The pronounced increases 
in NLO response on moving along the series 2^4 are con- 
sistent with their steadily extending it-conjugated structures. 

In order to determine the most appropriate theoretical 
method for treating Ru ammine complexes, we have tried 
various approaches to model the MLCT excited states. Pre- 
vious gas-phase DFT calculations using the B3P86 function- 
al with the LANL2DZ basis set proved qualitatively useful, 
but lack quantitative accuracy.' 21 ' A similar basis set combi- 
nation LANL2DZ(Ru)/6-31G* or 6-31 + G*(H,C,N,0) has 
been used by Inerbaev et al. ,' 23 ' and in a recent investigation 
by Zhang and Champagne.' 24 ' Here, we use a larger basis set 
(LANL2TZ(f)) for Ru in DFT calculations, as well as differ- 
ent larger basis sets for various high-level ab initio methods. 
We provide a more detailed account of the methodology in 
the Computational Methods. 

Comparing the new computational results with the experi- 
mental data shows that TD-DFT with the hybrid functionals 
B3LYP, B3P86 or M06 yields very good values for the first 
dipole-allowed excited state (FDAES) of 3, when MeCN 
solvent is included via the polarisable continuum model 
(PCM; Table 1). Selected simulated spectra are shown in 
Figure 2, together with the experimental spectra of com- 
plexes 3 and 4 as their PF 6 salts. Comparisons with data 
published previously' 23,24 ' show that although using a larger 
basis set for Ru does not affect the qualitative picture of the 
excitation properties, some significant quantitative differen- 
ces are observed. 

The excitation into the FDAES consists for 3 and 4 nearly 
exclusively of the HOMO^LUMO transition, whereas for 
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Table 1. Selected data calculated for complexes 2-4 by using DFT and ab initio meth- 
ods,'" 1 together with previously reported measured data. 



Method 


^max 

[nm] 


fc 


l max 
[nm] 


fc 


X max 
[nm] 


/» 


B3LYP 


439 


0.28 


449 


0.24 


439 


0.86 


djLY r/solvent 1 1 


456 


0.49 


578 


0.35 


641 

40gl c l 


0.69 
1.32' c ' 


M06/solvent 


472 


0.48 


599 


0.37 


659 


0.71 


B3P86 


437 


0.28 


454 


0.25 


444 


0.81 


rjJroO/SOlVclll 


4jj 


U.3U 


JO / 


U.JO 


4U/ 1 1 


Xi.IXi 

1 otM 


M06L/solvent 


455 


0.46 


788 


0.39 


937 


0.64 


CAM-B3LYP 


404 


0.27 


332 


0.17 


356 


1.81 


CAJVL-B JL.Y r/SOlvent 


434 


0.48 


379 


0.42 


416 


1.68 


Wriy /AiJ 


408 


0.27 


333 


0.12 


353 


1.82 


LC-coPBli 


374 


0.25 


277 


0.42 


329 


1.20 


LC-a)PBE/solvent 


393 


0.37 


305 


0.42 


347 


1.74 


RI-CC2/def2-TZVPP 


615 


0.42 


421 


0.34 


389 
389 


0.53 
0.57 


RASPT2(18,18) [d| 






437 


0.22 






RASPT2(18,18)/MeCN [d] 






585 


0.21 






B3P86/LANL2DZ [e) 






482 


0.37 


473 
366 


0.79 
1.06 


B3P86/LANL2DZ/MeCN |f| 






608 




693 




B3P86/6-31G*/LANL2DZ(Ru)/ 






546 


0.29 


607 


0.57 


MeCN [8] 










386 


1.12 


measured' 11 ' 


540 




590 


0.20 


584 
354 


0.43 
0.58 



[a] Molecules 2 and 3 were optimised with B3LYP/6-31G**/LANL2TZ(f) (Ru), 
whereas 6-31G* was used for 4. All optimisations were performed in the gas phase, 
except for 2, which was optimised in the gas phase and in H 2 Q [b] MeCN used as sol- 
vent for 3 and 4, H,0 used for 2. [c] For the second dipole-allowed excited state 
(SDAES). [d] State-averaged RASPT2. [e] Data taken from ref. [21]. [f] Data taken 
from ref. [23]. [g] Data taken from ref. [24]. [h] For [2][C10 4 ] 3 in H,0 (a very weak 
NIR band at l raax = 855 nm is observed also), |17e| [3][PF 6 ]„ |191 and [4][PF 6 ] 3 ; [21] A m „ 
values in MeCN at room temperature; directly corresponding /„ values are unavaila- 
ble, so those quoted are in PrCN at 77 K (only slight variations due to changing the 
solvent and temperature are expected). |20 ' 



2 it is a HOMO-l^LUMO transition. The orbitals, shown 
in Figure 3 for B3LYP as an example, demonstrate clearly 
that in each case a charge transfer from Ru to the organic 
ligand is involved, confirming the expected MLCT charac- 
ter. To better quantify these transitions, we have applied 
a model developed by Bahers et al. recently for the semi- 
quantitative analysis of photoinduced CT processes, which is 
based on the differences in electron density between the 
ground and excited states.' 25 ' Figure 4 shows the spatially-re- 
solved density differences upon excitation for the three com- 
plexes. This model also affords estimates of the transferred 
charge g cx and the CT distance r cx , which are 0.288 e/ 
1.571 A for 2, 0.953 e/4.218 A for 3, and 0.978 e/5.799 A for 
4, evaluated for B3LYP. To put these data into perspective, 
we note that the corresponding values for the prototypical 
CT molecule 4-nitroaniline in MeCN are 0.62 e/2.72 A. |25c ' 
The calculated r CJ values agree relatively well with the ef- 
fective (localised) electron-transfer distances r ab calculated 
from A ol « ab /e, where A 0L « ab is the dipole -moment change be- 
tween the diabatic states involved in the MLCT transition. 
Based on Stark spectroscopic measurements, respective r ab 



values of 3.6 and 5.7 A are determined for [3][PF 6 ] 3 
and [4][PF 6 ] 3 in PrCN glasses at 77 K.' 20 ' 21 ' 

The predicted l max values agree less well with 
those measured for 4, and the observed blue shift in 
the MLCT band on moving from 3 to 4 is not re- 
produced in solution. However, this blue shift is 
predicted by the gas-phase results, albeit with a de- 
crease in the overall accuracy of the A max values, 
reminiscent of previous gas-phase B3P86/LanL2DZ 
calculations.' 21 ' Our new calculations using B3LYP, 
B3P86 or M06 give rather large discrepancies be- 
tween theory and experiment for 2. Nonetheless, it 
is gratifying that the results obtained for 3 by using 
these functionals are very similar to those from 
high-level restricted active-space SCF second-order 
perturbation theory (RASPT2) with a very large 
basis set and large active space (see the Computa- 
tional Methods for details), both in solution and the 
gas-phase. The long-range corrected functionals 
(LRCFs) CAM-B3LYP,' 26 I LC-coPBE' 27 ' and 
wB97XB' 28 ' give A max predictions less accurate than 
those obtained when using B3LYP, B3P86 or M06 
for all three complexes (see also Figure 2 b). These 
results are quite surprising because one of the rea- 
sons LRCFs were introduced was specifically to im- 
prove descriptions of charge-transfer excitations by 
standard DFT functionals. ' 26 ~ 2S| However, our results 
concur with other recent studies on organotransi- 
tion metal compounds.' 29 ' For example, Escudero 
and Gonzalez' 290 ' found poor performance by 
CAM-B3LYP and LC-mPBE for MLCT excitations 
in fra«s-[RuCl 2 (2,2'-bpy)(CO) 2 ], when compared to 
experimental and RASPT2 results, with more accu- 
rate predictions from hybrid functionals. The rea- 
sons for such unsatisfactory performance of LRCFs 
are unclear at present. 

Although only gas-phase results are available for the fully 
ab initio and reasonably high-level resolution of identity ap- 
proximate coupled-cluster (RI-CC2) method, they differ 
substantially from those derived from DFT and RASPT2, 
and notably overestimate 2 max when compared with experi- 
ment for 2. For 4, RI-CC2 predicts two very close-lying tran- 
sitions (A max = 389 nm) of comparable intensity, which is also 
not predicted by any other method. The relative failure of 
RI-CC2 may derive from the rather large values of the 
diagnostic for the CC2 GS wavefunction (0.317, 0.124 and 
0.115 for 2, 3 and 4, respectively). Values of D t above 0.050 
may indicate a large multi-configurational character of the 
GS, for which the single-reference method CC2 is not suita- 
ble.' 30 ' However, without taking into account the apparently 
large solvent effect, it is difficult to judge the performance 
of RI-CC2 with confidence. As a general point, we note that 
DFT takes into account the non-dynamical correlation, al- 
though only partially and in a non-systematic way.' 31 ' This 
aspect may explain why the method performs better for 
complex 3 than for 2, which according to the D x diagnostic 
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Figure 2. a) Normalised electronic absorption spectra of 2 (simulated) 
and 3 (simulated and experimental) in solution; all the simulated spectra 
are convolutions of the computed B3LYP values with a Gaussian func- 
tion with cr=2000 cm -1 , b) Experimental and simulated electronic ab- 
sorption spectra of 4; both CAM-B3LYP and B3LYP simulated spectra 
are shown. The experimental data were obtained with the complex salts 
[3][PFJ, or [4][PF 6 ] 3 in MeCN. [w ' 211 



from CC2 may have stronger multi-configurational charac- 
ter. This point is still under investigation. 

Further tests using different optimised geometries show 
that the geometry does not have a large influence on A max or 
/ os . This observation validates the comparisons between the 
results of the RASPT2 calculations for 3, which for compu- 
tational efficiency reasons used a geometry of C s symmetry, 
and the data obtained from the other methods that used 
a Cj symmetry. The relative geometry independence holds 
also for the NLO properties of 3 and 4, but not for those of 
2. Therefore, all of the properties of 2 were computed with 
the geometry optimised in the corresponding environment 
(gas-phase or water solvent). 

Considering that the NLO properties are very dependent 
on a good description of the excited-state manifold, it seems 
that, apart from highly accurate (and computationally ex- 
pensive) multi-configurational methods, only the properties 
calculated with B3LYP, B3P86 or M06 may be at least ap- 
proximately reliable. Therefore, parameters for the GS and 




HOMO-1 



LUMO 




HOMO 



LUMO 



HOMO 



LUMO 



Figure 3. Orbitals (B3LYP) involved in the main transition from the GS 
to the FDAES for 2 (top), 3 (middle) and 4 (bottom). 




Figure 4. Plots showing difference electron density between the GS and 
the FDAES of 2 (top), 3 (middle) and 4 (bottom); green denotes positive 
differences, whereas negative differences are in red; isocontour values 
±0.001. 



FDAES derived by using these functionals for the com- 
plexes embedded in a solvent continuum are collected in 
Table 2, together with CAM-B3LYP results. All of the (hy- 
per)polarisabilities were computed by finite-field derivatives 
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Table 2. Further data calculated for complexes 2-4 by using DFT and RASPT2 methods,' 8 ' together with previously reported measured data. 



Complex 


Method 


State 








Yzzzz 




\>l a zz 












[au] 


[au] 


[10 3 au] 


[10 4 au] 


[au] 


[au] 


[au] 


[10 3 au] 


2 


B3LYP 


GS 


0.37 


336 


6 


-4 






2.71 


3.8 






FDAES 


1.23 


224 


17 


22 


0.86 


-112 








B3P86 


GS 


0.42 


337 


6 


_5 






2.74 


3.7 






FDAES 


1.25 


229 


17 


11 


0.83 


-108 








M06 


GS 


0.57 


346 


5 


-15 






2.74 


2.6 






FDAES 


1.11 


246 


17 


3 


0.54 


-100 








CAM-B3LYP 


GS 


-0.16 


310 


9 


39 






2.62 


8.0 






FDAES 


1.97 


139 


9 


279 


2.13 


-171 








measured' c ' 


GS 










K0 81 |d| 


(2201 






3 


B3LYP 


GS 


-0.69 


519 


-39 


133 






2.57 


-48 






FDAES 


-8.21 


238 


220 


-4900 


-7.52 


-281 










Relaxed S,' e ' 


-5.85 


295 


566 








1.43 






B3P86 


GS 


-0.76 


529 


-42 


138 






2.63 


-50 






FDAES 


-8.09 


205 


221 


-4700 


-7.33 


-324 








M06 


GS 


-0.79 


542 


-44 


1425 






2.70 


-52 






FDAES 


-7.69 


320 


233 


— 1500 


-6.90 


-222 








CAM-B3LYP 


GS 


-0.23 


405 


_11 


5000 






2.29 


-16 






FDAES 


-7.59 


834 


136 




-7.36 


429 








RASPT2(18,18) 












-6.44 




2.03 


-26 




literature' 1 ' 


GS 










-6.16 




2.54 


-33 




measured' 8 ' 


GS 










(5 4V d| 


(2301 


2.10 


28 


4 


B3LYP 


GS 


-2.35 


1035 


-116 


6800 






3.83/4.21 


-186/-214 






FDAES 


-13.06 


1158 


917 


-21000 


-10.71/-3.26 


123 








B3P86 


GS 


-2.42 


1057 


-126 


7300 






3.87 


-232 






FDAES 


-12.94 


1025 


912 


-23720 


-12.52 


-32 








M06 


GS 


-2.45 


1063 


-127 


7240 






3.94 


-195 






FDAES 


-12.50 


1265 


809 


-28900 


-10.05 


198 








CAM-B3LYP 


GS 


-1.69 


766 


-19 


1354 






4.79/1.80 


-70/-76 






FDAES 


-7.74 


3200 


-51 




-6.05/-5.42 


2434 








literature' 1 ' 


GS 










-9.54 




3.58/3.78 


-130/- 82 




measured' 8 ' 


GS 










(8.8)' 111 




3.10/2.70 


126 



[a] All data calculated in H 2 0 (2) or MeCN (3 and 4); 3 and 4 optimised in vacuum, 2 optimised in H 2 0. [b] 2SA = two-state approximation; 3SA = 
three-state approximation, p calculated from 6S" 1 A oi Uj(// oj ) 2 /(£' oi ) 2 where fi & is the transition dipole -moment and E tti is the transition energy from the GS 
to the FDAES (i = l) or SDAES (i = 2); n=l (2) for p 1SA (/8 3s a); an additional term for 3SA containing pi u was neglected because this property computed 
in the gas-phase is almost zero. The value for 4 is the total obtained by applying the 2SA to the two low-energy absorption bands separately (the lowest 
energy band with MLCT character involves the FDAES, whereas the other band has intraligand charge-transfer character and involves the SDAES).' 20 ' 

[c] Data for [2][BF 4 ] 3 taken from ref. [18]; numbers in brackets indicate properties measured in (or derived from data measured in) PrCN glasses. 

[d] Only the magnitude was determined, [e] Optimised structure in S! state, [f] Data taken from ref. [24]. [g] Data for [3/4][PF 6 ] 3 taken from ref. [20] and 
[21]. 



of the TD-DFT excited-state dipole -moments, and are thus 
static (zero-frequency) values. 

Most of the FDAES properties were computed for the op- 
timised GS geometries.' 32 ' Nevertheless, in order to assess 
the possible effects of excited-state structural relaxation on 
the properties, an excited-state optimisation of 3 in the gas- 
phase was carried out, and properties were calculated also 
for the resulting geometry (see the Computational Methods 
for details). 

The CAM-B3LYP results are again clearly very different 
from those obtained with the other functionals. Previously 
measured A 0L a z and A 01 a- Z values were obtained by Stark 
spectroscopy at 77 K in glassy PrCN solution, and are thus 
not directly comparable to computed values in liquid 
MeCN. Nevertheless, the latter are of the right order of 
magnitude when using B3LYP, B3P86 or M06. Interestingly, 
a zz is predicted to be smaller for the FDAES than for the 
GS for 2 and 3, but the reverse is apparent for 4 when using 
B3LYP or M06. [33] The calculations with B3LYP, B3P86 or 
M06 all predict the observed increases in A 01 ju z on moving 



along the series 2^4, and the Stark-based value of 8.8 au 
for 4' 21 ' is closest to that obtained when using M06. 

The GS first hyperpolarisability of complex 3 in the two- 
state approximation (j3 2SA ) from RASPT2(18,18) is very 
close to the value deduced from experimental HRS data 
(28472 au).' 19 ' The corresponding values from the DFT 
methods are considerably larger; the reason for this can be 
traced back primarily to a substantially larger transition 
dipole-moment predicted by the DFT methods, which is 
squared in the equation for /? 2 sa> an d is approximately 1.2- 
1.3-times larger than the RASPT2 result. Also, the DFT cal- 
culations give slightly larger dipole-moment changes than 
does RASPT2. 

We find that the presence of diffuse functions in the basis 
set used to describe the first and second row atoms, which is 
generally of great importance for the NLO properties of the 
GS, has little influence on the NLO properties of the 
FDAES. The same can be said, to a lesser degree, of polar- 
ised functions, which generally are considered to be impor- 
tant for reliable descriptions of excited states. These obser- 
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vations appear quite unusual. On the other hand, the basis 
set LAN2TZ(f) for Ru is about the minimum necessary for 
a reliable description of both the GS and FDAES proper- 
ties. This last factor explains the main differences between 
our data and those obtained by Inerbaev et al. and by 
Zhang and Champagne, who used the smaller basis set 
LAN2DZ.' 23 ' 24 ' 

When considering the responses calculated by using 
B3LYP, B3P86 or M06, these increase substantially (as ex- 
pected) on increasing the Jt-conjugation path-length along 
the series 2^4, for both the GS and FDAES species 
(Table 2). Also, substantial enhancements are found for the 
FDAES with respect to the GS in each complex, with 
changes of similar magnitude predicted when using the 
three different functionals. However, the /3 values calculated 
by using CAM-B3LYP show less consistent trends. Hence, 
although the GS Ru n -containing chromophores possess 
large NLO responses, the activity is even larger for the 
MLCT excited states that formally contain Ru m coordinated 
to a reduced pyridyl ligand radical. The predicted excited- 
state enhancement of /J becomes more significant as the jt- 
conjugated system extends, being about threefold for 2, ap- 
proximately fivefold for 3 and about sevenfold for 4. For 3, 
geometry relaxation within the FDAES (the S t state) leads 
to further considerable enhancement of /J, whereas a increas- 
es slightly. 

The cubic NLO properties, denoted by the second hyper- 
polarisability y, are generally more difficult to compute reli- 
ably with our method, and thus should be considered as 
more approximate. Even so, they are also generally in- 
creased in the FDAES (Table 2). However, unlike for /J, the 
results vary between the functionals B3LYP, B3P86 and 
M06. Using B3LYP or B3P86 yields enhancements of y in 
each case, but these are much larger for 3 (ca. 34-37-fold) 
when compared with 2 or 4 (ca. two- to sixfold). In contrast, 
the calculations with M06 predict an excitation-induced in- 
crease in y for 4 only (ca. fourfold), whereas no change is 
evident for 3 and a decrease is found for 2. Unfortunately, 
no experimentally measured y values are available for 2-4 
to allow comparisons with theory, because these complexes 
are of interest primarily for their quadratic NLO properties. 
The general result that both /3 and y are larger in the 
FDAES is consistent with previous studies involving other 
types of chromophore.' 8 " 11 ' 

Conclusion 

Compared with previous computational studies on the elec- 
tronic excitation and GS quadratic NLO properties of Ru" 
ammine complexes, we have used a larger basis set 
(LANL2TZ(f) for Ru) in DFT calculations. Also, we have 
used different larger basis sets for various high-level ab 
initio methods. TD-DFT calculations with the hybrid func- 
tionals B3LYP, B3P86 or M06 and a MeCN PCM yield very 
good agreement with the experimental spectrum for excita- 
tions to the FDAES of 3. The same methods give a lower 



degree of matching with the experimental data for the short- 
er or longer complexes, underestimating l max for 2, but over- 
estimating it for 4. In each case, the calculations confirm the 
expected MLCT character of the FDAES. A model devel- 
oped by Bahers et al. affords electron-transfer distances that 
agree well with those determined via Stark spectroscopic 
measurements on [3][PF 6 ] 3 and [4][PF 6 ] 3 previously. Notably, 
the results obtained for 3 by using B3LYP, B3P86 or M06 
are very similar to those from RASPT2 with a very large 
basis set and large active space, both in solution and the 
gas-phase. Surprisingly, A max values predicted by the long- 
range corrected functionals CAM-B3LYP, LC-coPBE and 
wB97XB are less accurate than those obtained with B3LYP, 
B3P86 or M06 for all three complexes. Gas-phase results 
from the fully ab initio RI-CC2 method differ substantially 
from those derived from DFT and RASPT2, possibly due to 
the relatively high multi-configurational character of the 
GS. 

Considering both the GS and FDAES, polarisabilities and 
hyperpolarisabilities were computed by finite-field deriva- 
tives of the TD-DFT excited-state dipole-moments, by using 
the functionals B3LYP, B3P86 or M06. The A 01 /t z and A 01 a zz 
values are of magnitude similar to those measured previous- 
ly by Stark spectroscopy in the frozen solution state at 77 K. 
The GS /3 2SA value for 3 from RASPT2(18,18) is very close 
to that derived experimentally by HRS, whereas the corre- 
sponding DFT-based values are considerably larger, primari- 
ly due to substantially larger predicted fi n values. Somewhat 
surprisingly, the presence of diffuse and polarised functions 
in the basis set used to describe the first and second row 
atoms has little influence on the NLO properties of the 
FDAES. As expected, the /? responses calculated by using 
B3LYP, B3P86 or M06 increase markedly as the Jt-conjuga- 
tion extends along the series 2^4, for both the GS and 
FDAES species. Also, substantial enhancements are found 
for the FDAES with respect to the GS in each complex, 
with the three different functionals predicting changes of 
similar magnitude. The excited-state enhancement of /? in- 
creases as the n-conjugation extends, from approximately 
threefold for 2 to about sevenfold for 4. Although more ap- 
proximate and of less interest for the complexes studied, the 
computed y values also generally increase in the FDAES, 
but the results vary between the functionals B3LYP, B3P86 
and M06. 

In summary, state-of-the-art theoretical methods show 
that MLCT excitation of Ru 11 ammine chromophores leads 
to large increases in molecular NLO responses. Therefore, 
such complexes hold promise not only as redox-switchable 
species, but also for ultrafast photoswitching of bulk NLO 
effects in appropriate organised media, such as LB thin 
films. 



Computational Methods 

The molecules 2 and 3 were optimised with the 6-31G** basis set for 
CHN and the LANL2TZ(£) ECP/basis set for Ru at the DFT level with 
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the B3LYP functional. For 4, the 6-31G basis was used for H. Molecule 2 
was optimised in the gas-phase and in aqueous solution, by using the 
PCM, and the respective structures used for corresponding phase compu- 
tations; for 3 and 4, gas-phase optimised structures were used for all cal- 
culations except for the RASPT2 computations (see below). 
Excited-state optimisation of 3 in the gas-phase was carried out with 
three different basis sets for the atoms H, C, N, O: 6-31G**, 6-31 + G** 
and 6-31 + +G**, while LANL2TZ(f) was applied to Ru throughout. In 
order to be comparable with the other calculations, the properties were 
computed with the 6-31G** basis set and the PCM/MeCN model. The re- 
sults can be grouped into two different sets according to the basis sets un- 
derlying the excited-state optimisation. With the 6-31 + G**-optimised 
structure, a small Stokes shift, A st , of 0.28 eV was obtained, with a large 
oscillator strength for the transition to So (/ Os = 0.43). This is approximate- 
ly in line with the results reported by Zhang and Champagne,' 241 (A st = 
0.25 eV, / os = 0.31), who apparently also used this basis set for the excited- 
state optimisation, although with LANL2DZ for Ru, the B3P86 function- 
al and the PCM/MeCN model already applied during the optimisation. 
However, contrary to their result of the relaxed state being S b we find it 
to be S 3 , which is also the Franck-Condon state. Usually, such a highly 
excited state would relax into Si by internal conversion before geometry 
relaxation is completed. For the 6-3IG** and 6-31 + +G**-optimised 
structures, on the other hand, the relaxed state becomes Si, showing 
large A St values (1.2 and f.l eV for the 6-31G** and 6-31+ +G** opti- 
mised structure, respectively) and small /„ values for dipolar transitions 
into So (0.05 for both structures), which is in reasonable agreement with 
the absence of luminescence observed experimentally. Consequently, we 
used the 6-31 G**-optimised structure to compute the excited-state prop- 
erties. However, it should be noted that the optimisations of the second 
group were not fully successful in terms of the required criteria applied 
in Gaussian 09, even after about one hundred cycles. The best structures 
were finally used, where three of the four criteria were fulfilled. 
TD-DFT with several functionals was employed to compute transition 
moments and excited-state dipole moments, which in turn were used for 
finite-field numerical differentiations to obtain the hyperpolarisabilities 
in the FDAES. The dipole moments of the charged species are reported 
with respect to the centre of nuclear charge. The GS properties were cal- 
culated analytically or by finite-field differences from GS dipole mo- 
ments. Only the dominant long-axis component (oriented in z-direction) 
was computed in all cases. Most of the PCM excited-state calculations 
were performed by using the state-specific description for the non-equi- 
librium solvation during a vertical excitation,' 34 ' although the differences 
when compared with the linear response non-equilibrium approxima- 
tion' 35 ' were generally negligible. The transition energies were measured 
from the electronic GS minimum, and no vibrational averages were 
taken into account. It should be noted that the comparison of computed 
vertical excitation energies with the experimental absorption maxima is 
only an approximation. A more accurate method would take the Franck- 
Condon factors (overlap between the vibrational wavefunctions of 
ground and excited electronic states) into account also. However, such an 
approach is computationally much more expensive, 1 ' 6 ' and has thus not 
been used here. The useable field strengths for the finite-field numerical 
differentiations had to be adapted for each molecule and functional indi- 
vidually: the lowest useful field strengths depended on the magnitude of 
the properties, while the high-field limit was set by field-induced state 
mixing and switching. Gaussian 09' 37 ' was used for these calculations. 
The RI-CC2 calculations were carried out with Turbomole.' 38 ' For the 
RASPT2 calculations on complex 3, its structure was first optimised at 
the PBE0/TZVP/MWB-28(Ru) level in C s symmetry with Turbomole. Al- 
though this structure is not a minimum, the symmetry was required to 
make the multi-configurational calculations manageable. Tests with this 
structure at the DFT level showed that the influence of the non-mini- 
mum structure on the properties of interest is minimal. State-averaged 
RASPT2 calculations with extended ANO-rcc basis sets, 
([8s7p6d3f2glh] (Ru), [4s3p2dlf] (C,N), [3s2pld] (H)) were then per- 
formed on this structure with Molcas.'^' Extensive tests were performed 
to find the optimal active space. These led finally to the following parti- 
tion : the RAS2 space contains eight orbitals, five of Ru 4d parentage, 



two of their o bonding counterparts and the Jt-orbital involved in the 
FDAES. RAS1 contains four Jt-orbitals and RAS3 three jt* and three 
Ru 5d orbitals. Single and double excitations were allowed out of RAS1 
and into RAS3, while all kinds of excitations were allowed in the RAS2 
space. 

For the three-state approximation of the first hyperpolarisability of 4 
shown in Table 2, the excited-state dipole-moment between the FDAES 
and the second dipole-allowed excited state was computed with the 
Dalton2011 package, in the gas phase. As the resulting value was very 
low, the corresponding term in the sum-over-states expression was ne- 
glected. Finally, the difference electronic density calculations were per- 
formed by using the CT program.' 41 ' 
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